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O ■ The Bethe ansatz equations of the 1-D 577(3) Hubbard model are systemat- 

ically derived by diagonalizing the inhomogeneous transfer matrix of the XXX 

qs ' model. We first derive the scattering matrix of the SU(3) Hubbard model 

through the coordinate Bethe ansatz method. Then, with the help quantum 
inverse scattering method we solve the nested transfer matrix and give the eigen- 
values, the eigenvectors and the Bethe ansatz equations. Finally, we obtain the 
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1 Introduction 

Strongly correlated electron systems have been an important research subject in the 
condensed matter physics and mathematical physics. One of the significant models 
is the 1-D Hubbard model for which the exact solution was first given by Lieb and 
Wu[0. Based on the Bethe ansatz equations (Lieb-Wu's equations), the excited spec- 
trum (spin excitation and charge excitation) was discussed in Refs.0, ^, [|, [5], ^|. The 
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magnetic properties of the 1-D Hubbard model at zero temperature were investigated 
in Refs.[[7|, ||, §], and the thermodynamics of the model was studied in Refs. |T0|, p] , |T2|| 
with the help of string hypothesis. The critical exponents of the system was found in 
Refs. fl3|, (L4]]. All physical properties can be discussed in the framework of Bethe ansatz 
equations. Although there are lots of works on Hubbard model, the integrability of 



the 1-D Hubbard model was finished until 1986 by Shastry|15|, Olmedilla and Wadati 
jlBl . Moreover, the eigenvalue of the transfer matrix related to the Hubbard model was 
suggested in Ref.[15] and proved through different method in Ref. [17, 18]. Besides, the 
integrability and the exact solution of the 1-D Hubbard model with open boundary 



condition have been investigated by several authors. [19, 20, 21 



Recently, based on the Lie algebra knowledge, Maassarani and Mathieu constructed 
the Hamiltonian of the SU(n) XX model and shown its integrability p2|. Consider- 
ing two coupled SU(n) XX models, Maassarani succeeded in generalizing Shastry's 
method to construct the SU(n) Hubbard model [23|. Furthermore, he found the re- 



lated .R-matrix which ensures the integrability of the one- dimensional SU(n) Hubbard 
model|M]].(It was also proved by Martins for n = 3, 4.[25|, and by Yue amd Sasaki [25| 
for general n in terms of Lax-pair formalism.) However, the eigenvalue and the eigen- 
vectors of the SU(n) Hubbard model have not been given yet. 

In this paper, we first apply the coordinate Bethe ansatz method to construct the 
wave function and the scattering matrix of the £77(3) Hubbard model in section 2, then 
with the help of the Yang-Baxter relation we list the algebra consisting the elements 
of the transfer matrix in section 3, which is the key relations for finding the solution of 
the nested transfer matrix. In section 4 we apply the algebraic Bethe ansatz method 
to discuss the eigenvalue and the Bethe ansatz equations of the model. Then the exact 
solution of the SU(3) Hubbard model was given out. Under the thermodynamical 
limit, the explicitly analytic form of the ground state energy is given in section 5. In 
section 6, we make some conclusions and list some questions to be considered. 

2 The Coordinate Bethe Ansatz 

The Hamiltonian of the SU(n) Hubbard model is 

L n-1 tj 2 L 

tt \ "■ \ V rpna rpan , rpan rpna , rpna rpan , rpan rpna , ^ \ "■ ft r~i l-\\ 

H ~ 1^ l^\ h a,k h a,k+l + Ar.kAr.fc+l + K,kK,k+l + K,kK,k+l + —r~ 2^ ^a,k^r,k (1J 



k=la=l 



fc=l 



where U is the coulumb coupling constant, and E" k (a = a,r) is a matrix with zeros 
everywhere except for an one at the intersection of row a and column (3: 



(E^) lm = 5?5i 



(2) 



The subscripts a and k stand for two different E operators at site k(k = 1, • • • , L). The 
n x n diagonal matrix C is defined by C = J2 a <n E aa — E nn . We have also assumed 
the periodic boundary condition E%+ L = E% . 

It was shown that the Hamiltonian (P has a (stt(n — l)©tt(l)) CT ©(sw(n— 1)©m(1)) t 
symmetry. The generators are 

n 
T a/3 _ J^ R a/3 

k=l 

and 



Jf = E E t ( 3 ) 



K a = ^2C a>k , a,p=l,---,n-l, a = a,r. (4) 



fc=i 



It is worthy to point out that the system enjoys the 5*0(4) symmetry when n = 2 [28]. 
But the generator are different from Eq. (|J). In this paper, we just limit our attention 
on the SU(3) model. 

Before proceeding the coordinate Bethe ansatz approach, we first begin by intro- 
ducing some notations of our analysis. In SU(3) Hubbard model, there are two types 
of particle named a and r, and each type of particles can occupy one of two possible 
states. We denote |0) the vacuum state, |1) and |2) the two particle states respectively. 
Under the appropriate basis 



2) = 1 , |0) = (5) 






we can prove that E a3 and E 3a act as a creating and a destroying operators of \a) 
state respectively. 

In the coordinate Bethe ansatz method, the eigenstates of the Hamiltonian can be 
assumed as 

k*o> = E fzrz:° e:? x1 e:? X2 . . . #$*„ io>, (6) 



xi<X2<---<x^ =l 



where 



ai a 2 ---a No _ y^ , a Q\ a Q2 '" a QN (i 7 i \ 

J(T1<T2-<TN — 2-^ e P e Q Aa Ql a Q2~ a QNr ) \ Kp Ql^ P Q 2 ' " ' ^PQnJ 

P,Q 

No 

x6(x Ql <x Q2 <--- <XQ No )exp(i^2k Pj x Q .), (7) 

i=i 

Here ckj = 1, 2(i = 1, 2, • • • , N ) stands for the different particle states, Xi the position of 
the particle, Cj = cr,T the type of ith particle. The summation P and Q are taken over 



all permutations of iVo momenta kj and iVo coordinates Xj respectively. The symbols ep 
and eQ are the parities of two kinds of permutations. Substituting the wave function 
into the Schrodinger equation 

H\1> No ) = E\1> No ), (8) 

we can get 

a::: ; ^';;;(- • • , h, %,•••) = s%» (sin k h sin fc^iCg&cG ■ ■ , %, h, ■ ■ ■). (9) 

with Sa-ai (smki, sin kj) being the two-particle scattering matrix: 

S3* (sin k„ sin t,) = ™*.-^+^ (10) 

1 J sin /Cj — sin kj + 27 

where 7 = ^, P<^ J is the direct product of two kinds of permutation operators and 
permutes the particle styles and particle states simultaneously. The energy of the 
Hamiltonian on this wave function is 

i^ 7 

E = 2j2co8ki + i-(L-2N ). (11) 

i=i ^ 

For the convenience, we denote by Co the amplitude A... |"j ' CT /,'- (•■-,%, fo, ■ • •). When 
the j'-th particle moves across the else particle, it gets an ^-matrix <%((& — g/). Here 
qj = sin kj. The periodic boundary condition leads the following constrains on the 
amplitude Co 

Sjj+iSjj+2 ' ' ' SjN SjiSj 2 ■ ■ ■ Sjj-iCo = e 1 J Co- (12) 

Eq.fll^D is similar with the Yang's eigenvalue problem[^7|. Its solution will give out the 
Bethe ansatz equation. 

3 The foundamental commutation rules 

In the above section, we have obtained the scattering matrix. It can be proved that 



they satisfy the Yang-Baxter equation p7 



Sji(qj, qi)S jk (qj, qk)Si k (qi, q k ) = Si k (qi, q k )S jk (qj, q k )Sji(qj, qi). (13) 

Define C = (Sjj+iSjj+2 • • • «S , jjv )~ 1 Co> we can rewrite the eigenvalue problem (|P2|) as 

S jl S j2 ---S m C = X J C = e(q j )C (14) 



We introduce an auxiliary space r and S T j = S T j(q — qj), and define the monodromy 
matrix 

T T (q) = S t1 S t2 ■ ■ ■ S tNo (15) 

Obviously, S r j(qj) = P T j and 

tr T T T ( qj )=X r (16) 

Since SV« satisfies the Yang-Baxter equation, we can easily prove that the monodromy 
matrix T T (q) also satisfies the Yang-Baxter relation 

S T r'(q, q')T TJ (q, qj)T T/j (q', %) = T T ,j(q', qj)T Tj (q, qj)S TT ,(q, q'). (17) 

From the point of view of a vertex model, we can interpret the matrix S T j as the 
vertex operator, the matrix S TT > as the /2-matrix. So the i?-matrix is an 16 x 16 matrix 



Rji(lj,Qi) 



where 
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(19) 



(20) 



where dj, tj are two different kinds of Pauli matrixes on j-th space. 

In addition to the L-operator and the i?-matrix the existence of a local reference 
state is another important object in the quantum inverse scattering program. In our 
case, we can define the local reference state by 



The action of L-operator on this state has the following property: 



/ |0)i 





C TJ (q)\0), 



Uh 







|0>i 
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q + %r I 



(21) 



The global reference state |0) is then defined by the tensor product of local refer- 
ence states, i.e. |0) = ITJLi ®\fyj- This state is an eigenstate of the transfer matrix 
which is the trace of monodromy matrix. In order to construct the other eigenstates, 
it is necessary to seek for an appropriate representation of the monodromy matrix. By 
this we mean a structure which is able to distinguish creation and annihilation oper- 
ators as well as possible hidden symmetries. The property of L-operator suggests the 



monodromy matrix to take the following form[17 



Tr{q) 



B(q) B(q) F(q) 
C(q) A(q) B*(g) 
C(q) C*(g) D(q) , 



m\ 



where B(g),C*(g) and B*(g),C(g) are two component vectors with dimensions 1x2 
and 2x1 respectively. The operator A (q) is a 2 x 2 matrix and we shall denote its 
elements by A a b(q). The remaining operators B(q),C(q), F(q), D(q) are scalars. 

In the framework of the above partition the eigenvalue problem for the transfer 
matrix becomes 

[B(q) + J2 A aa (q) + D(q)]\$) = A(g)|*>. (23) 

a=l 

where A(q) and |$) correspond to the eigenvalue and the eigenvector respectively. From 



Eq. fl2ip , we know the action of the monodromy matrix on the reference state: 

S(g)|0) = |0), D(q)\0) = u(q)\0), 

A aa (q)\0)=u(q)\0), a = 1,2. 

C(g)|0) = 0, C*(g)|0> = 0, 

C(g)|0) = 0, A ab (q)\0) = 0, a^b. 

where the function u>(q) is defined by 



(24) 



A'o 



LO 



1~1j 



(?) = n 

J=i 1~1j + H 



(25) 



The operator B(g) and F(q) play the role of creation operators over the reference state 

|o). 

To make further progress we have to recast the Yang-Baxter algebra in the form of 
commutation relations for the creation and annihilation operators. We shall start our 
discussion by the commutation rule between the operators B(g) and B(p): 



B(q) ® B(p) = [B(p) ® B(g)] ■ f(g,p) 
where r(q,p) is an auxiliary 4x4 matrix given by 

/ 1 \ 

a(q,p) b(q,p) 

b(q,p) a(q,p) 

\ X ) 

and the functions a(q,p) and b(q,p) are defined by 



(26) 



f(q,p) 



(27) 



a(g,p) 



Q=3(g,P) 



6(g,p) 



QJ2(g,j?) 
ai(g,p) 



(28) 



We can see that this auxiliary matrix r(q,p) is precisely the rational i?-matrix of 
the isotropic 6-vertex model or the XXX spin chain. It is very easy to see out that 



a(q,p) = 1 -b(q,p), 



(29) 



To solve the eigenvalue problem (24) we need several other commutation rules, 
especially the commutation rules between the diagonal and the creation operators. 
The commutation relations between the diagonal and creation operator B(g) are 

A(q) B(p) = ^4t B (P) ® A (^)l • ? (*>P) ~ ^t4 b (q) ® A(p), (30) 

B(g)B(p) = 2^B(p)B(g) - ^BpM, (31) 

«2 ® 2 (p,q) 

D(q)B(p) = B(p)D(q) + ?^M).F(p)C*(q) - ?^±F(q)C*(p) (32) 

a 2 (q,p) a 2 {q,p) 

We also need the following commutation rule 

C*(g) <g> B(p) = ^^[B(p) ® CT(g)] • f (g,p) - ^4 B (?) ® C*(p) (33) 

0:2(9, p) a 2 (g,p) 

We have set up the basic tools to construct the eigenvectors of the eigenvalue 
problem ( fffi ). In the next section we will show how this problem can be solved with 
the help of the commutations rules (|26|), d27|) , (|30|) - C|33|) . 



4 The eigenvectors and eigenvalue construction 

The purpose of this section is to solve the eigenvalue problem for the transfer matrix. 
We shall begin by considering the construction of an ansatz for the corresponding 
eigenvectors. 

4.1 The eigenvalue problem 

The eigenvectors of the transfer matrix are in principle built up in terms of a linear 
combination of the products of many creation operators acting on the reference state, 
which are characterized by a set of rapidities parameterizing the creation operators. 
First, we define the eigenvector for an arbitrary A^-particle state: 

l$*i({Pl})> = *J\ri(Pi>P2, • • • ,p Nl ) ■ ^|0) (34) 

where the mathematical structure of the vector ^ Nl (p 1 ,p 2 , ■ ■ ■ ,PjvJ will be described 
in terms of the creation operators. We denote the components of vector T by jF ai " ajv i 
which will be determined later on, where the index a, runs over two possible values 
a* = 1,2. 



To construct the eigenvectors of the transfer matrix, it is sufficient to look for 
combinations between the operators B(g) and F(q). The experience on constructing a 
few particle states suggests us the eigenvector of iVi-particle to be 

*jvi(Pi, • • -,Pni) = B (fi) ® *jvi-i(P2, ■ • -iPm)- (35) 

Here we have formally identified $o with the unity vector. The detailed construction for 
less than 3 particles was given in Appendix B. This vector has the following symmetry: 

$jvi(Pi, • • -,Pj-i,Pj, ■ ■ ■ ,PNt) = *JVi(pi, • • -,Pj,Pj-i, ■ ■ • ,PNi_) ■ rj-ij(Pj-i,Pj), (36) 

where the subscripts in rj-ij(j)j-.i,pj) emphasize the positions in the iVi-particle space 
V\ ® • • • <8> Vj-i <g) Vj <8> • • • <g) V^ on which this matrix acts non-trivially. Here we have 
already assumed that the (N\ — l)-particle state was already symmetrized. 

Applying the diagonal elements of monodromy matrix on this eigenstate, we can 
get 

tJi^yPiA) 
I>(g)|«M{ft})) = w(g)|$ Wl ({p,})). (38) 

E4(g)fe(W)> = ^(Q)U^r^^ 1) (<i,{pimN 1 ({pi})) 

£i a2(g,Pi) fc= i,t^ a 2(Pi,Pfc) 

x|*g;_ 1 (g,p„fe})}. (39) 

here 

l*$-i(g,Pi, fe})) = B(g) ® ^-iCpi, • • • ,fc, • • • ,p Nl )dP(pi, {pi}) • -F|0). (40) 
and 

4-1 

Oj (ft, {p*}) = II h,k+i(j>k,Pi) (41) 

fe=i 

where the symbol pj means that the rapidity pi is absent from the set {pi, • • • ,Pn ± }- 

The terms proportional to the eigenvector |$zvi({p*})) are denominated the wanted 
terms because they contribute directly to the eigenvalue. From Eqs. ([37])- ([fl]), we can 



9 



directly get the eigenvalue of iVi-particle state 

a(«, {«}) = fi ^#4 + «(«) + «(«) n ^') (4 fe ». 



= i«2fe,g) 



= i«2(g,Pi 



(42) 



The remaining ones are called unwanted terms and can be eliminated by imposing 
further restrictions on the rapidities pi. These restrictions, known as the Bethe Ansatz 
equations, are 






(43) 



In fact, the undeterminated eigenvalue A^\q, {pi}) must satisfy the following aux- 
iliary problem 



il-'-ajv-. 



(44) 



where the inhomogeneous transfer matrix T^(q, {pi}) is 

T W (q, {R})?;^ = CX(<l,PiK b d l(q,P2) ■ ■■rt 1 1 cl bNl (q,PN 1 ). (45) 

This result is the direct extensions of that obtained from the two-particle state. 
The Bethe ansatz equations and the eigenvalue still depend on this additional auxiliary 
eigenvalue problem. 

4.2 The eigenvalues and the nested Bethe ansatz 

The purpose of this subsection will be the diagonalization of the auxiliary transfer 
matrix T^'(q, {pi})- First, we write the transfer matrix T^\q, {pi}) as the trace of the 
following monodromy matrix: 



T {1 \q, {pi}) = ^(^(^paO^waw^'^i- 1 ) ' ' •^iui^Pi)- 

-a) 



(46) 



where A^ is the auxiliary space. The L-operator Cj^ (q,Pi) is related to the auxiliary 



matrix f(q : Pi) by a permutation operator 

/ 1 \ 

b(q,pi) a(q,pi) 



£%i(Q,Pi 



a(q,pi) b(q,pi) 
V 1 ) 



(47) 
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Comparing with 6- vertex model, one can find that this problem is equivalent to the 
6- vertex model. So we can take the results of Ref.[jl7H directly. The eigenvalue of the 
auxiliary problem PJ3I is: 



N 2 1 JVi N 2 1 

A (1) (g, {Pi}{s k }) = II r ( r + II Hq,pi) II r, r- ( 48 ) 

k=1 b{s k ,q) l=1 fc=i%,Sfc) 

and the parameters {sk} satisfy the following restrictions: 

Nl - N2 b(si. s) 

n^*,ft)=n?r 5if T. k=i,...,M. (49) 

1=1 j=l 0(Sj, Sk) 

Now we use the auxiliary eigenvalue expression to rewrite our previous results for 
the eigenvalues and the Bethe ansatz equations. Substituting the expression ( fUjj) in 
Eqs. (02,03) and using the relation (E5I) we obtain that the eigenvalue is 



w r XJ U ft fa-g + ^7) , -S (g-ft) , S (g-^j 

x(n ( v p;+ ; 7) n (g r g+ ; 7) + n (9 7 3fc+ : 7) k5o) 

[jii (?-pj) fc =i (a*-?) fe =i (?-<*) ' 

In terms of the above eigenvalue, the periodic boundary condition (|l|) implies 

e^ L = Afo, to}, k}) = n [ sin ^:f ( 5l ) 

£* (sm % - P; + zf ) 
The second set of Bethe ansatz equations Eq. (|43|) for the rapidities pz changes into 

i=ii& (Pi -Pi- l V j= i {pi - sin fcj - if) fc= i (pi - «fc - *i) 



(52) 



and the third set of Bethe ansatz equations Eq. (PS) for the rapidities s& is 



n frr?^ n i 4 ' ft ; i7 ' (*=i,2,-,a«. (5 3) 

Here we have used the shifted parameters p~; = pi — z| adn s~fe = s k — ij to bring our 
equations to more symmetric forms. 
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5 The ground state 

In the proceeding section, we have obtained the final Bethe ansatz equations of the 
SU(3) Hubbard model. In this section, we will use them to analyse the ground state 
property of the model. 

By taking the logarithm of the Bethe ansatz equations (|5l| ) - (|53|) , we can get the 
following three sets of equations (In this section, we assume U > 0): 

kjL = 2-kIj + 2 f] tan" 1 f ^hzA \ , (54) 

1=1 \ 2 / 

arJl 4- 2 t tan- (*Z*) = 2 £ tan"' (ELZ^l) + 2 £ tan"' ( «Z*) , 

i=l,*^ V 7 / j=l \ 2 / fc=l V 2 / 

(55) 
2^tan- 1 f^^)=2 £ tan" 1 ( *LZ*i\ + 2 nP kl (56) 

1=1 \ 2 / i=l,»^A: \ 7 / 

where Z,, J;, P^ are integers or half-odd integers. Under the thermodynamic limits the 
Bethe ansatz equations for the ground state change into 

v ; 7-s 7 2 + (A - A') 2 

, y - dk+ f E 2 47 ;' A '> 2 ^, (58) 

-Q7 2 + 4(A-sinA;) 2 7-b 7 2 + 4(A - A') 2 

2?rT (A) + /* _W) dA , = f W) 

v ; 7-B 7 2 + (A - A') 2 7-B 7 2 + 4(A - A') 2 

where Q,B and .E are determined by the conditions 

j_ Q p(k)dk = ^, (60) 

B iVi 

<r(AWA = — , (61) 

-b L 

J_ E r(A')dA' = ^. (62) 

12 



The function p(k),a(A),r(A') are the distribution functions of real parameters kj,pi 
and s~fc respectively. 

Eqs. (0)"O) nave a unique solution which is positive for all allowed Q, B and 
E. jf, (i = 0, 1, 2) is a monotonically function of Q, B and E respectively. Thus the 
ground state is characterized by Q = n, B = E = oo. 

After taking Fourier transforms of the above equations we can obtain the result of 
the distribution functions 

1 cosk f o° l + e -7H _ sinfc 



p(k) = / ' - ■ v e- wsink Mu)duj, (63) 

PV ; 2vr 2vr J-oo 1 + 2cosh( 7 M) 0V ; ' V ; 

1 /■«> 2cosh(||w|) 



27T J-oa 4 COsh (f |W|) — 1 



r(A') = ^-/°° , , 2 ^, h - JbMe-^W (65) 

27T j-oo 4 cosh (^|u;|) — 1 

From Eqs.(61)-(66) we have A^q = -^, this means that all lattice site is filled by one 
particles, N\ = |iV , and N 2 = ^N . 

The ground state energy then is given by 

/OO ry 

p(k) cos kdk - -L 
-oo Z 

roc I _|_ e ~7|w| ~ 

= -2L/ .— 3 JoMJxM^-lL. (66) 

J-oo u[l + 2cosh(7|o;|)J 2 

6 Conclusions 

The main purpose of this paper was to investigate the eigenvalue of the SU(3) Hub- 
bard model. We have succeeded in finding the eigenvalue of the Hamiltonian ([I]) and 
obtained three sets of Bethe ansatz equations. Based on the Bethe ansatz equations, 
we have found the explicit expression of the energy and the distribution functions of 
the rapidities corresponding to the ground state for positive U. 

One important question is to study the excitation spectrum and the low tempareture 
thermodynamics of the model for both positive and negative U. The Bethe ansatz 
equations given in present paper will play an key role. As we know, the negative U has 
a distinguished properties from positive U case. The solution structure of Bethe ansatz 
equations corresponding to the ground state and excited spectrum for attractive case 
(U < 0) are different from those for repulsive case. It is worthy to be studied in future. 



In Ref . ||24j| , the L-operator and i?-matrix were given, therefore, one can define the 



transfer matrix as done in usual Hubbard model ( SU(2) case). But how to find 

13 



the eigenvalue of the transfer matrix is unknown. It may be solved by using the 



method proposed in Ref.||17||. We will consider this problem late. Another natural 
generalization the present paper is to find the eigenvalue of the Hamiltonian ([!]) for 
general n, it will be considered elsewhere. 

7 Appendix A: 

First the one-particle eigenstate can be assumed to be 

W = ECWO). (A.i) 

x=l 

Applying the Hamiltonian ([!]) into this ansatz, we can obtain 

9U 
EQ(x) = Q(x + 1) + Q(x - 1) + —{L - 2)Q(x), 1< x < L. (A.2) 

The solution of above equation is 

f%(*) = K\^ kx ~ K\-^ ik \ (A.3) 

with the eigenenergy 

9U 
E = 2cosfc+— (L-2), (A.4) 

The periodic boundary condition governs the momentum to be 

kj = ^, (i = l,.-.,L-l). (A.5) 

For the two-particle state, the eigenstate is assumed as 

life> = E Q?M>^)K? xl K?jo)- (A.6) 

xi,a;2=l 

From the Schrodinger equation (|j) we have 

9U 
Ef%£(x llXa ) = niTM + 1, **) + CS'fo - 1. **) + — ( L ~ 4)/ CT T 2 2 (^i,x 2 ) 

+f%%(xi,X2 + 1) + C^l, *2 - 1), Xl ^ X 2 , 

9U 
EQ^{x u x 1 ) = f%£(x 1 + l,x 1 ) + f%%(x 1 -l,x 1 ) + - r f%g(x 1 ,x 1 ) 

+f££(xi,xi + 1) + f^(x 1 ,x 1 - 1), Xl = x 2 . (A.7) 

14 



The fermionic property of particle requires the antisymmetry of the wave function. 
Therefore, we can assume the wave function f^" 2 (xi,x 2 ) to be 

{Aa 1 a 2 (l, h, \„ik 1 x 1 +ik 2 x 2 _ Aa 1 a 2 ( jL h, \„ik 2 x 1 +ik 1 x 2 „ ^ _ 

-Klal{ k i: k 2 )e iklX2+ik2Xl + A^(k 2 , ki)j k2X2 + iklXx , x x > x 2 y } 



Substituting the above ansatz into the Eq. ( |A.7| ), we can obtain the eigenvalue of the 
Hamiltonian ([!]) 

9U 

E = 2(cos fci + cos k 2 ) + —(L- 4). (A.9) 

Imposing the continuous condition of the wave function 

^lim + /££(*!, x 2 ) = hm_ f£%(x 1: x 2 ), (A.10) 

we have 

A^lihM) = S%%(sink u sink 2 )A%%(k 2 ,h), (A.ll) 

where 

QQl a 2 1 ■ , • . x sin fci - sin fc 2 + H^V? , A -, ^ 

^<x lCT2 ( sm fc i' sm fc 2j = :— : :— : — , . (A.12) 

sm fci — sm k 2 + 27 

is the scattering matrix of two particles with labels (<ti, «i) and (02, 02), and 7 = -^-, 
P"^ 2 is the direct product of two kinds of permutation operators. 

8 Appendix B: 

Following the procedure of Ref . [I7j , the one-particle eigenstate is 



I^i(Pi)> = Bfo) ■ ^|0) = 5 a ( Pl )^|0). (B.l) 

where the repeated index means sum. 

With the help of commutations (Bn)-(B2) and the properties (|2"E|), (P7|), we find that 
the one-particle state satisfies the following relations 

BMI-MP,)) = ^4*^0) " ^B(,) . m , (B.2) 

i>(<z)|$i(pi)) = w(g)|$i(pi)), (B.3) 

a=l «2W,PlJ 

-^O^^B^.^IO). (B.4) 
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The terms proportional to the eigenvector |3>i(pi)) denominate the eigenvalue 

Afa, *) = ^4 + «(«) + w (,)afe£iA«(,, *). (B.5) 

«2(Pi,g) "2(9, Pi) 

where A^(g,pi) is defined by the one-particle auxiliary eigenvalue problem 

T {1 \q^) a c l b a\F bl = r a c l h al(q,Pi)F hl = A (1) (?,Pi)^ cl (B.6) 

The remaining ones will be eliminated by imposing the Bethe Ansatz equation 

w(pi) = I- (B.7) 

The two-particle state should be a composition of two single hole excitations and 
a local hole pair. The former is made by tensoring two creating operators B(g) while 
the latter represented by F(q). Thus, we construct the two-particle vector as: 

$ 2 (p!,p 2 ) = B(pO <g> B(p 2 ) +^(Pi)5(p 2 )^ 0) (Pi,P 2 ), (B.8) 

where #q (pi,p 2 ) is an arbitrary function to be determined. The vector £ plays the 
role of an "exclusion" principle, forbidding two same particles at the same site. 

Putting the diagonal elements of transfer matrix on this state, we find that the 
function g$ (pi,p 2 ) must be zero such that the two-particle state is the eigenstate of 
the transfer matrix. In other words, using the definition ([34] ) we have 



|$ 2 (p l5 p 2 )) = BifrjBjfa^lO). (B.9) 

The diagonal operators acting on this state gives the following relations: 

-E^H n ^^l*i 1, (?,P.;te})).(B.io) 

~i Oi2\JPi, q) kJi^i a 2 {Pk,Pi) 
D(q)\^ 2 (p 1 ,p 2 )) = w(?)|$ 2 (pi,P2)>, (B.ll) 

J2A aa (q)\UPi,P2)) = ^(g)f[^7^A( 1 )(g,{p,})|$ 2 (p 1 ,p 2 )) 

o=i i=i«2ig,Pij 

^ .a 3 (g,pi) i ai(pi,p fc ) (1) , f u 

<=i «2(g,p 4 ) fc= i,t^«2 (Pi,Pfc) 

x|*?W,«)> 5 (B.12) 
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where we have used the two-particle auxiliary relation 

T {1 \q, {pi})H%F a ™ = rg£(g,ft)Off>Pa)^ ai = ^^ {ft})** 9 *, (B.13) 

and 

\y ( i\q,Pi, {Pi})) = B(g) <8> ^(p^O^fe, fe}) • JF|0>, k ^ i, (B.14) 

The vanishing of unwanted terms gives us the Bethe ansatz equation 

1 = ^ ^^) A a) fe , {M ), j = i >2 (B.15) 



Finally, From Eqs. (|B . 1 0|) - ( p3 . 1 2|) we can directly read the eigenvalue of two-particle 
state 

a(«, to}) = n ^4 + *(«) + <"(«) n ^r^AWfe, { W }). (B.i6) 

Before constructing the general multi-particle state, let us first analyse the symme- 
try of the two-particle state. From the above analysis, it is easy to prove the two-particle 
vector satisfying the following exchange property 

*2(pi,P2) = *2(P2,Pi) • r(p 1 ,p 2 )- (B.17) 

In principle, such symmetrization mechanism can be implemented to any multi-particle 
state. This indeed help us to handle the problem of constructing a general multi-particle 
state ansatz. We construct the three-particle state as 

*s(pi,P2,P3) = B(pi) ® B(pa) ® B(p3), (B.18) 

satisfying 

*3(P1,P2,P3) = *s(P1,P3,P2) • ^2,3(P2,P3)- (B.19) 

We also can rewrite the three-particle state in terms of the following recurrence 
relation: 

*3(Pi 5 P2,Ps) = B(pi) ® $ 2 (p 2 ,P 3 )- (B.20) 

This expression is rather illuminating , because it suggests that a general iVi-particle 
state can be expressed in terms of the N\ — 1-particle state via a recurrence relation. 
This recursive way not only help us to better simplify the wanted terms but also makes 
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it possible to gather the unwanted terms in rather closed forms. The action of diagonal 
operators on the three-particle state is 

B(q)\$ 3 ( Pl ,p 2 ,p 3 )) = n^|«,(p 1)ft) p 5 )) 

y^ « 3 fa,g) ^ ai(pk,Pi) 
^i<X2(pi,q) k= [ l k ^a 2 (pk,Pi) 

x|*?W,{pi})>, (B.21) 

D{q)\$a(pi,P2,P3)) = u(q)\® 3 (p u p 2 ,p 3 )), (B.22) 

J2A aa (q)\®s(pi,P2,P3)) = u(q)i[^^\AU(q,{pi}m3(pi,P2,P3)) 

V* f \ a s(,Q,Pi) TT a l(Pi'Pk) A (l)/ r U 

xl^^Pi,^})}. (B.23) 

In general, the knowledge of the 2-particle and the 3-particle results suggests the 
behaviour of the A^-particle state. By using the mathematical induction we are able 
to determine the general structure for the multi-particle states to be Eq. ([34]) 
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